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Context. The cosmic-ray ionization rate fcR of dense molecular clouds is a key parameter for their dynamics and chemistry. 
Aims. Variations of <fcR are well established, but it is unclear if these are related to source column density or to Galactic location. 
Methods. Using the APEX telescope, we have mapped the 364 GHz line of H3O* in the Sgr B2 region and observed the 307 GHz line at 
selected positions. With the IRAM 30-m telescope we have observed the H'**0 203 GHz line at the same positions. 

Results. Strong H3O* emission is detected over a ~3x2pc region, indicating an II3O* column density of lO'^-lO" cm"- in an 18" beam. The 
HaO^ abundance of ~3xlO"' and H^O^/HtO ratio of ~l/50 in the Sgr B2 envelope are consistent with models with <J'cr~4x10"'^ s"', 3x lower 
than derived from H^ observations toward Sgr A, but lOx that of local dense clouds. 

Conclusions. The ionization rates of interstellar clouds thus seem to be to first order determined by the ambient cosmic-ray flux, while propa- 
gation effects cause a factor of ~3 decrease from diffuse to dense clouds. 
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1. Introduction 

The ionization rate of interstellar clouds i s a key parameter for 
, their dyn amics and their c hemistry (see ICaselli & WalmslevI 
■ I2OOII and lVan der Takl2006l for reviews). In dense star-forming 
regions, the bulk of the ionization is due to cosmic rays, at a rate 
of fcR~3xlO s"' as derived from HCO^ and CO emission 
line observations of the envelopes of young massive stars up to 
a few kpc from the Sun (Van der Tak & van Dishoeck 2000). 
This result is in good agreement with measurements of the local 
flux of low- energy co smic rays with the Voyager and Pioneer 
spacecrafts llWebbei< fl998). However, mm-wave observations 
of DCO+ and HCO+ toward the nearby da rk cloud L DN 1544 
indicate a decrease of ^cr by a factor ~10 dCaseniet al. 2002). 
while observations of absorption toward the P er diffuse 
cloud indicate an increase of similar magnitude ( iLe Petit et alJ 



However, in the case of Sgr B2, using the DCO^/HCO^ ratio 
is complicated by the uncertain D abundance, while comparing 
HCO^ wi th CO is difficult as the molecules may not trace the 
same gas (! jaca et all 199^. which calls for other methods. 

The high proton affinity of water makes hydronium (HsO^) 
a key ion in the oxygen chemistry of dense clouds. Its sub- 
millimeter rotation-inversion transitions may be used to mea- 
sure ^CR and trace H2O and O2 w hich are unobservable from 
the ground. Pionee ring work by IWootten et al.l (Il99ll) and 
'Philli ps et aP (Il992!) has demonstrated this potential with ob- 
servations of strong H30^ emission toward Sgr B2 (OH)'. 
Unfortunately, the atmospheric transmission from Mauna Kea 
is rarely good enough to observe HsO^, so mapping is not fea- 
sible. 

The Chajnantor site greatly facilitates observation of the 
H3O+ 7/f=3+-22 hne at 364.7974 GHz, which has an up- 



12004 ). Even larg er enha ncements are found for the Galactic per level energy of 139 K and an Einstein A-coefficient of 
Center, where Ok a et alJ lEo05) derive ^cr~10"'^ s"' from H+ 2.8x10-"^ s"'. The Jk^I'-I^^ line (v=307. 1924 GHz, £„=80K, 



observations toward Sgr A. An enhanced ionization rate in the 
inner 250 pc of our Galaxy is expected from its strong X-ray 
and radio synchrotron emission, but this value is at the upper 
limit based on the observed temperatures of the Sgr B2 clouds 
jGiisten et alJH985.) . It is unclear if these variations in ^cr are 
due to propagation effects (absorption or scattering of cosmic 
rays) or to variations in the cosmic-ray flux with location in 
the Galaxy. To resolve this issue, an estimate of the ioniza- 
tion rate of dense gas in the Galactic Center is urgently needed. 



A=3.5xlO" s" ) is a probe of denser gas, because the emis- 
sion competes with fast pure inversion decay, unlike for the 
364 GHz line. This paper presents new observations of these 
lines toward the Sgr B2 region. Combined with observa- 



' This position between the Sgr B2 (M) and (S) cores was used for 
early line surveys of the Galactic Center fe.g.. lCummins et'al]ll986h 
but doe s not corre spond to any object, not even the centroid of the OH 
masers JGaume & ClausserLl990-) . 
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Table 1. Measured line parameters at selected positions with 
Icr errors in units of the last decimal in brackets. Upper limits 
are Icr limits on in K on 1.0 km s"' channels. 



Position 


Qr(J2000) 5(J2000) Jr,^bdV 


Vlsr 


AV 




hh mm ss " ' " Kkm s"' 


km s"' 


km s"' 



H3O+ 364 GHz: 



SgrB2(M) 17:47:20.2 -28:23:05 42.0(13) +62.0(6) 17.3(16) 

SgrB2 (OH) 17:47:20.8 -28:23:32 8.7(4) +62.7(3) 15.3(8) 
GCS 3-2 17:46:14.9 -28:49:43 <0.09 
GCIRS3 17:45:39.6-29:00:24 <0.10 





H3O+ 307 GHz: 


SgrB2(M) 


27.1(52) +64.6(14) 16.6(21) 


SgrB2(OH) 


7.7( 7) +64.9(13) 24.0(22) 




Hf 203 GHz: 


SgrB2(M) 


20.5(3) +63.3(1) 13.2(1) 


SgrB2(OH) 


<0.03 



tions of the O 3i3-2o2 Une at 203.4075 GHz (£„=204 K, 
A=4.9xl0"^ s"'), the data are used to estimate ^cr in the 
Sgr B2 region. Due to the high critical densities of the lines, 
our observations are not sensitive to the e xtended low-density 
envelope s een in ab sorption lines of water ( Co mito et alJ2003l) 
and H30^ (iGoicoechea & Cernicharo.2001.) . 

2. Observations 

In July 2005, we observed the 364 GHz hue of H3O+ towards 
selected sources in the Sgr B2 and Sgr A clouds using the 
APEX telescope (Glisten et al., this volume)^. The front end 
was the facility APEX-2a receiver built at Onsala (Risacher et 
al., this volume); the back end was an FFT spectrometer built 
at the MPIfR providing 8192 channels over a bandwidth of 
1.0 GHz (Klein et al., this volume). A map of Sgr B2 at 10" 
spacing (Fig. Q was made by integrating 1 minute per point, 
which for system temperatures of 340 - 360 K gave a noise 
level of 0.2 K per 0.8 km s ' channel. Deeper integrations (2- 
5 min) were made at the (N), (M) and (OH) positions. The line 
was also observed toward two positions in the Sgr A cloud, 
CGS 3-2 and GC IRS 3, where stror ig mid-infrar ed absorp- 
tion suggests a high ionization rate tokaetalJ2 0()?). The non- 
detections at these positions probably reflect a lack of dense 
gas along these lines of sight. Follow-up observations of the 
307 GHz line of H^O^ were performed in November 2005. 
Since this line lies very close to a CH3OH line, several local 
oscillator settings were used to verify the sideband origin of 
the detected features. In this case, the system temperature of 
300 K resulted in an rms noise level of 90 mK on 0.24 km s"' 
channels after 10.4 minutes of on-source integration. 

All data were obtained using position switching with a 
throw of 3' in Az. Test spectra taken at the Sgr B2 (OH) off'set 
position indicate that the emission is <10% of the strength at 
the target position. The beam size of APEX is ~18" at 364 GHz 
and ~21" at 307 GHz. Data at both frequencies were calibrated 

^ This paper is based on data acquired with the Atacama Pathfinder 
Experiment, which is a collaboration between the Max-Planck- 
Institut fiir Radioastronomie, the European Southern Observatory, and 
the Onsala Space Observatory. 
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Fig. 1. Distribution of velocity-integrated H30^ 364 GHz emis- 
sion toward Sgr B2, observed with the APEX telescope. Filled 
dots indicate detections; open dots denote upper limits. The 
greyscale runs from 1.9 to 35.5 in steps of 3.7 Kkm s"'. 

onto Tmb scale using forward and main beam efficiencies of 
97% and 74%. From a comparison of calibrations using hot 
and cold loads and using a water vapour radiometer, we es- 
timate the calibration uncertainty to be wlO%. For reasonably 
strong lines (T^\, I lOxrms), as is the case here, spectral rms has 
a negligible contribution to the intensity uncertainty. 

Observations of the Hj^^O 203 GHz line were performed 
with the IRAM 30-m telescope^ in February 2006. The facility 
receiver A230 was used as front end and the VESPA autocorre- 
lator as backend. At this wavelength, the telescope has a beam 
size of 12" and a main beam efliciency of 57%. System tem- 
peratures were 500-600 K and integration times 10-30 min, 
resulting in rms noise levels of ~0.1 K on 0.20 km s ' chan- 
nels. The line was detected at the Sgr B2 (M) but not at the 
(OH) position. Calibration wa s verified to 10% on Sgr B2 (N) 
which was observed before bv lGensheimer et al] ( Il996l) . 

3. Results 

Emission in the H^O^ 364 GHz line is detected over 
an (80x40)" region in Sgr B2 (Fig. corresponding to 
(3.3x1.6) pc at a distance of 8.5 kpc. The overall distribution 
of H30^ 364 GHz in the Sgr B2 region follows that of many 
other molecular lines when mapped at comparable reso lution 
(e.g.. 'Mar tin-Pintado et al1ll990Hde Vic ente et al."l997^. The 
distribution of H30'^ also follows that of the (sub-)millimeter 

^ IRAM is an international institute for millimeter-wave astron- 
omy, co-funded by the Centre National de la Recherche Scientifique 
(France), the Max-Planck-Gesellschaft (Germany) and the Instituto 
Geografico National (Spain) 
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Frequency (GHz) 

364,80 364,78 364,76 




Table 2. Column densities of H2O, H3O+, and H2 in the Sgr B2 
region. 



60 _ 80 

V,,sH (km s ') 
Frequency (GHz) 

307,22 307,20 



SgrB2(M) 



SgrH2(0H) (x2) 
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Fig. 2. Spectra of the H3O+ 364 GHz Une (top) and 307 GHz 
Une (bottom) toward the Sgr B2 (M) and Sgr B2 (OH) posi- 
tions. 



dust continuum emission d Lis et al.ll99lh . a common tracer of 
the H2 distribution. In contrast, the H2O 183.31 GHz emission 
is strongly peaked at the co res, dropping to ; 10% of its peak 
value at the (OH) position jCernicharo et alJl2006l) . This be- 
haviour does not necessarily reflect the H2O column density 
distribution because of the maser nature of the 183 GHz line. 

Table ^ lists the results of fitting Gaussian profiles to the 
spectra after subtracting a polynomial baseline. The 307 GHz 
line strength toward OH may be overestimated due to blending: 
fixing Ay to 15-16 km s"' reduces the line flux by 25 -30%. 
The st rength of this line agrees with the data of Phillips 'et alJ 
(Il992h . while the 364 GHz line is a factor of 2 weaker. The 
steep emission gradient (Fig.Q makes the measurement sensi- 
tive to pointing offsets. In the APEX data, H^O^ and other lines 
peak at the M core, validating the pointing to <5". 

The 307 and 364 GHz spectra show many lines in ad- 
dition to H3O+. While the Sgr B2 (OH) spectra show <10 
mostly weak (r,nb!0.5 K) lines over the 1 GHz bandwidth, 
the Sgr B2 (M) spectrum shows >20 features with Tmb ! 1 K. 
The spectra taken in the vicinity of Sgr B2 (N) are so confused 
that the H30^ line parameters cannot be reliably extracted. 
Identification and interpretation of these other lines will be 
done elsewhere in the framework of a spectral line survey pro- 
gram (Belloche et al., in prep.). 

4. Column densities of HsO^ and I-I2O 

Figure |2 shows spectra of the two H30^ lines at the posi- 
tions where 307 GHz observations have been made. At both 
positions, the intensities of the two H30^ lines are compa- 
rable. Since the lines are optically thin for any chemically 
reasonable H30^ abundance, the Une ratio of ~unity indi- 



Component 




A'(H30+) 


A'(H2)" 




lO'^cm-- 


10'^ cm-2 


1024 


Core (M) 


35 


13.6 


4.4 


Envelope (OH) 


0.15* 


3.7 


1.1 


Beam size (") 


12 


18 


20 



" : From CSO 350 fjm data assuming Tj=40 K 
and Ky=0.07 crn^g^^ (Lis, priv. comm.) 
From ^emkh^^^l^ 022© assuming H2O/H2 



lo- 



cates a high excitation temperature, Tex ! 50 K. Statistical e qui- 
librium calculations using molecular data from Schoier et alJ 
12005)4 indicate that at the conditions in the M core (r=200K, 
«(H2)=10^ cm"-*), collisions can sustain Tex^lOOK, consis- 
tent with the observed H30^ line ratio. However, the density 
(10^ cm"^) at the OH position is too low to thermalize the ex- 
citation at the kinetic temperature of 60 K. The observed H30"^ 
Une ratio therefore indicates that the H30^ excitation is driven 
by far-infrared pumping by the strong dust continuum radiation 
in the Sgr B2 region, as found before bv .Fhillips et al.. 
Models with r^=40-80K reproduce the observed Une ratio. 
Far-infrared pumping is also expected to play a major role to- 
ward (M), due to its high infrared luminosity, and is therefore 
included in the model. 

For these excitation conditions, the observed line strengths 
imply A?(p-H3O+)=1.83(±0.27)xl0'^ cm^^ for Sgr B2 (OH) 
and 6.81(+1.01)xl0'^ cm ^ for Sgr B2 (M). The total H3O+ 
column densities are twice these values, as the ortho/para (o/p) 
ratio tends to unity at high temperatures (; 100 K). These es- 
timates have uncertainties of a factor of ~2 due to the uncer- 
tain excitation and the difference in estimated values from the 
two lines. Since the lines are optically thin and the source sizes 
cannot be constrained, the derived column densities are beam 
averages. For the Sgr A sources, assuming Ay=10km s"', 
Tex^lOOK and o/p-1, the observed limit on the line flux of 
wl Kkm s"' corresponds to A^(H30^) <6xl0'-' cm"^. 

Dividing the above column density estimates by A^(H2) 
values based on (sub-)millimeter dust continuum mapping 
(Table |2j provides estimates of the column-averaged abun- 
dance in the two sources: A^(H30+)/A^(H2) - 3.1x10"'' in 
Sgr B2 (M) and 3.4x10"^ in Sgr B2 (OH). These estimates of 
the H30^ colu mn density and abundance are higher than pre- 
vious values bv Phi llips et alJ (Il992l) which were derived from 
the H30^ 396 GHz line assuming a density of 10^ cm for the 
Sgr B2 e nvelope, a value which l ikelv is an order of magnitude 
too low (Goico echea et alJ2004l) . 

To estimate the H2O column density from the IRAM obser- 
vations, we adopt the same excitation conditions as for H30^. 
To convert p-H^^^O into H2O, we adopt o/p=3 as applicable 
for war m gas and '^O/'^O=2 50 as applicable for the Galactic 
Center dWilson & Roodll994 . For the Sgr B2 (M) position, we 
estimate A^(H2O)=3.5xl0'^cm"^, while the upper limit at the 
(OH) position impUes A?(H2O)<1.7xl0i'^cm-2 (Icr). Within 
the uncertainty of a factor of two, our estimate of A^(H20) at the 



|http : //www ■ strw . leidenuni v . nl/~moldata/radex . php| 
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M position is consistent with radiative tra nsfer modelUng of th e 
H^^^O 203 GHz Hne towards the N core dComito et alJl2003l) . 
The H2O column densities imply abundances of 1.2x10"^ and 
<1.7xl0"^ at the (M) and (OH) positions, respectively, consis- 
tent with the estimates of ~10~^ and ~10~^ frorn observations 
of the H2O 183 GHz line jCernicharo et alJ2006l) . 

5. Chemistry of H2O and H3O+ 

Table 12] summarizes measurements of the H2O and H30^ col- 
umn densities toward the Sgr B2 cores (represented by the M 
position) and the surrounding envelope (represented by the OH 
position). The H30^/H20 ratio is seen to be -1/50 in t he enve- 
lope, or about 20 times the value in local clouds dPhillips et alJ 
I1992 I). but drops by ~100 toward the cores. This drop is more 
than expected from enhanced recombination at higher densi- 
ties, and may be due to ongoing injection of H2O into the gas 
phase by evaporation of icy grain mantles, which also causes 
the rich (sub-)millimeter spectra of the (M) and (N) cores. 

To estimate the cosmic -ray ionization rate ^cr in the Sgr B2 
region from our observations, we have used the chemical 
model of Van der Tak & van Dishoeck ( 2000). Given estimates 
of temperature, density, ^cr, and the abundances of CO, O, N2 
and H2O, this model calculates the steady-state abundances of 
HCO^, N2H^ and HsO^, as well as the total electron fraction. 
The model neglects metals such as Mg, Fe and S as electron 
sources and PAHs as electron sinks; including these species 
may change the calculated electron fraction by 2-3 depend- 
ing on their abundance. The reaction rates in t he model have 
been updated to follow the UMIST99 values jLe T euff et all 
I2OOOI) . except for the dissociative recom bination r eactions of 
H+, HCO+ and N2H+ whe re th e values oflM cCall et al. (2003), 
Ide Boisanger et al.l (I199(tI) and lCieppert et aL (12004) are used. 
For the cores, we adopt T - 200 K and « = 10^ cm"^, 
whil e for the envelope, we u se T = 60 K and n - 10^ cm"^ 
(see IGoicoechea et alJl2004 . The CO abundan ce of 5x10'^ 
is from the SEST Hne survey (jNummelin et ai1 l2000l) . the O 
abundance of 1.5xlO~^ is based on the ISO-LWS line survey 
( IGoicoechea et al.ll2004l) . and we assume equal CO and N2 
abundances. 

The model reproduces the observed H30^/H20 ratio in the 
Sgr B2 envelope for ^cr~4x10"'^ s"'. The uncertainty is a fac- 
tor of 2 from the observed H30^/H20 ratio, and another fac- 
tor of 2 through the assumed CO, N2 and O abundances. For 
Sgr B2 (M), the model indicates a ~10x lower ^cr, which we 
consider a lower limit because the chemistry may not be in 
steady state. The time scal e to equilibrate th e H20/H30^ ra- 
tio is /ion ~ 30(«c) °^ yr dVan der Tak et alj Eb06) which for 
the model electron density of n^^O.Ol cm"-' equals 300 yr. 

The derived value of ^cr for Sgr B2 (OH) is lOx higher 
than the value for dense gas in the Solar neighbourhood (§ 1), 
but 3x lower thari indic ated by H^ observations of the Sgr A 
region l Okz. et al.l2005l) . We conclude that the ionization rates 
of dense molecular clouds are mainly determined by their lo- 
cation in the Galaxy through variations in the cosmic -ray flux 
by a factor of ~ 10. As a second order effect, ^cr is 3x lower in 
dense molecular clouds than in diffuse clouds, presumably due 
to cosmic-ray scattering dPadoan & Scalo 2005) . 



In the future, large-scale mapping of the Sagittarius molec- 
ular clouds in the H3O+ 364 GHz and H^^O 203 GHz lines, 
combined with constraints on the temperature and density 
structure of the clouds, will provide more detailed understand- 
ing of the nature of the interaction between cosmic rays and 
molecular gas. The IRAM and APEX telescopes are well suited 
to carry out such observations, especially if equipped with 
heterodyne array receivers. Observations of H2O and HsO^ 
ground-state lines with the Heterodyne Instrument for the Far- 
Infrared (HIFI) on the Herschel Space Observatory will allow 
us to extend these studies to regions of lower column density. 
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